clear
use simresults.dta

gen demmargin=districtdem-(100-10-districtdem)

gen demwinsym=1-repwinsym
gen demwinasym=1-repwinasym

collapse demwin*,by(demcandideo demmargin)
saveold simresults--unsmoothed.dta ,replace

* computing smoothed probabilities *
clear
use simresults--unsmoothed.dta
gen demwinsymsmooth=.
gen demwinasymsmooth=.
foreach pid of numlist -20(5)20 {
   lowess demwinsym demcandideo if demmargin==`pid',gen(a) nog 
   replace demwinsymsmooth=a if demmargin==`pid'
   drop a
   lowess demwinasym demcandideo if demmargin==`pid',gen(a) nog 
   replace demwinasymsmooth=a if demmargin==`pid'   
   drop a
}
replace demwinsymsmooth=1 if demwinsymsmooth>1 & demwinsymsmooth~=.
replace demwinsymsmooth=0 if demwinsymsmooth<0
replace demwinasymsmooth=1 if demwinasymsmooth>1 & demwinasymsmooth~=.
replace demwinasymsmooth=0 if demwinasymsmooth<0
saveold simresults--smoothed.dta ,replace
